arXiv:1504.03807v3 [cond-mat.mtrl-sci] 21 Apr 2015 


Dirac Line Nodes in Inversion Symmetric Crystals 


Youngkuk Kim/ Benjamin J. Wieder/ C. L. Kane/ and Andrew M. Rappe^ 

^ The Makineni Theoretical Laboratories, Department of Chemistry, 

University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA 
^Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6396, USA 

(Dated: April 22, 2015) 

We propose and characterize a new Z 2 class of topological semimetals with a vanishing spin-orbit 
interaction. The proposed topological semimetals are characterized by the presence of bnlk one¬ 
dimensional (ID) Dirac Line Nodes (DLNs) and two-dimensional (2D) nearly-flat surface states, 
protected by inversion and time-reversal symmetries. We develop the Z 2 invariants dictating the 
presence of DLNs based on parity eigenvalues at the parity-invariant points in reciprocal space. 
Moreover, using first-principles calculations, we predict DLNs to occur in CusN near the Fermi 
energy by doping non-magnetic transition metal atoms, such as Zn and Pd, with the 2D surface 
states emerging in the projected interior of the DLNs. This paper includes a brief discussion of 
the effects of spin-orbit interactions and symmetry-breaking as well as comments on experimental 
implications. 


A recent development in condensed matter physics has 
been the discovery of semimetallic features in electronic 
band structures protected by the interplay of symme¬ 
try and topology. A tremendous amount of progress has 
been made in materials with strong spin-orbit interac¬ 
tions, such as the surface states of topological insulators 
mm and topological crystalline insulators [3], as well as 
the gapless bulk states of Weyl and Dirac semimetals 01 - 
[6] . Related topological phenomena can occur in materials 
with vanishing (or weak) spin-orbit interactions [7]. In¬ 
deed, the prototypical topological semimetal is graphene 
[5] , which exhibits Dirac points that are robust to the ex¬ 
tent that the spin-orbit interaction in carbon is weak. In 
the absence of spin-orbit interactions, the Dirac points in 
graphene are topologically protected by the combination 
of inversion and time-reversal symmetries. 

In this paper we study a related phenomenon for three 
dimensional (3D) materials with weak spin-orbit inter¬ 
action. We show that the combination of inversion 
and time-reversal symmetries protects Dirac line nodes 
(DLNs), for which the conduction band and valence band 
meet along a line in momentum space, and we predict re¬ 
alistic materials in which they should occur. DLNs have 
been discussed previously in the context of models that 
have an additional chiral symmetry, which can arise on a 
bipartite lattice with only nearest neighbor hopping. In 
this case, the DLN can be constrained to occur at zero 
energy. However, chiral symmetry is never expected to 
be an exact symmetry of a band structure. We will show 
that despite the absence of chiral symmetry, the line node 
is protected, though it is not constrained to sit at a con¬ 
stant energy. We will show, however, that in the vicinity 
of a band inversion transition, a DLN can occur in the 
form of a small circle, whose energy is approximately flat. 
The presence of such a Dirac circle has interesting con¬ 
sequences for the surface states, and we show that on 
the projected interior of the Dirac circle, the surface ex¬ 
hibits a nearly flat band, which must be half-filled when 


the surface is electrically neutral. Such surface states 
could be an interesting platform for strong correlation 
physics. We introduce a class of materials and use first- 
principles density functional theory (DFT) calculations 
to show that they can be tuned through the band inver¬ 
sion transition and exhibit the predicted Dirac circle, as 
well as a more complex nodal structure. A similar DLN 
near an inversion transition has recently been predicted 
in a 3D graphene network [aiig. Recently, DLNs also 
have been proposed in systems with strong spin-orbit in¬ 
teractions, such as perovskite irridates mi III] and non- 
centrosymetric semimetals mi, but in those systems the 
mechanism of symmetry protection is different. 

We will begin by elucidating the topological constraints 
that inversion and time-reversal symmetries impose. We 
will then introduce a Z2 topological invariant (related 
to the invariant characterizing a 3D topological insu¬ 
lator) which signifies the presence of DLNs. We will 
then present DFT calculations on transition metal-doped 
CU3N that predict a Dirac circle, as well as nearly-flat 
boundary modes. We will then introduce a simple low- 
energy k ■ p model that explains the appearance of the 
Dirac circle at a band inversion, and allows for a simple 
description of the resulting boundary modes. 

Consider a 3D Bloch Hamiltonian H(k) that is invari¬ 
ant under inversion V and time-reversal T. In the ab¬ 
sence of spin-orbit interactions we may consider = 
-1-1. The occupied Bloch eigenstates are characterized 
by a Berry connection A(k) = — i ^^(u„(k)|VkU„(k)). 
V and T symmetries constrain the Berry phase, uj{C) = 
exp i A.- dk, on any closed loop C in momentum space, 
to satisfy w(C') = w(—C) and u){C) = lj(—C)*, respec¬ 
tively. It follows that loops C are characterized by a Z 2 
topological invariant uj(C) = ±1 [13]. The non-trivial 
loops uj(C) = —1 must enclose a degeneracy. In two di¬ 
mensions, this explains the symmetry protection of the 
Dirac points in graphene. In three dimensions, it guar¬ 
antees that a line of degeneracies must pierce any surface 
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bounded by C. 

The parity eigenvalues ^n(ra) = ±1 of the occupied 
Bloch states at the 8 parity-invariant momenta Tq pro¬ 
vide an important constraint that allows us to identify 
topologically protected line nodes. First, consider a time- 
reversal invariant loop Cab = Cab — Cab that connects Fa 
and Fb along two time-reversed paths Cab and Cab- In the 
supplementary material we prove that the Berry phase 
on this loop satisfies 


C^(Cab)=Ub; ea = n^"(ra)- (1) 

n 

If we now consider four parity-invariant points, the con¬ 
tour Cab — Ccd defines the boundary dSabcd of a surface 
Sabcd in momentum space. The Berry phase on dSabcd 
counts the number of DLNs N(Sabcd) that pierce that 
surface. We thus conclude that 

= UbUd- ( 2 ) 

Thus, when = ~I there must be an odd num¬ 

ber of DLN piercing any surface Sabcd, with the simplest 
case being just a single one. This relation is quite similar 
to the topological invariants (vq; i^it' 2 ^' 3 ) characterizing a 
(strong or weak) Z 2 topological insulator in the presence 
of spin-orbit interactions [El [E] ■ Indeed, in a topolog¬ 
ical insulator with = ~1 j when the spin-orbit 

interaction is turned off a DLN must appear, because the 
system can not be adiabatically connected to a trivial in¬ 
sulator. 

This connection to the parity eigenvalues suggests a 
route towards realizing the DLNs: Starting with a triv¬ 
ial insulator, invert a pair of opposite-parity bands. At 
the inversion transition a small Dirac circle will necessar¬ 
ily emerge and grow. In the following we will predict a 
class of real materials which exhibits this behavior, and 
analyze the low-energy structures which emerge. 

Searching for materials that consist of light elements 
and preserve T and V symmetries, we find that copper ni¬ 
tride, CU 3 N, a narrow-gap semiconductor (Eg ~ 0.3 eV) 
m, fosters DLNs near the Fermi level via an insulator- 
to-metal transition driven by doping transition metal 
atoms. Copper nitride, first synthesized in 1937 
is stable in air at room temperature with a cubic anti- 
ReOa structure in space group 221 (Pm3m). It contains 
a rather large void at the center of the cubic unit cell, 
as shown in Fig. [l] This void can host extrinsic atoms 
such as N [E] , Li [El HO] ) Pd PUHS] , Rh, Ru [E] , Zn, 
Ni, Cd [E], Cu [ElHO], Fe, Ti gT], Ag [55], La, Ce [E], 
as well as many other transition-metal atoms |30j . In 
particular Ni, Cu, Pd, An, Ag, and Cd m are found to 
drive an electronic transition in CU 3 N into a semimetal 
without breaking T symmetry |30j . by which we expect 
that DLNs form near the Fermi energy. 

To demonstrate the existence of DLNs in the transi¬ 
tion metal-doped CU 3 N, we perform first-principles cal- 



FIG. 1. Crystal structure of CusNX. X represents a transi¬ 
tion metal atom intercalated at the body-center of the cubic 
unit cell of CusN in an anti-ReOa structure. 


culations based on DFT. We employ the Perdew-Burke- 
Ernzerhof-type generalized gradient approximation |31j 
as implemented in the QUANTUM ESPRESSO pack- 
age |32j . Norm-conserving, optimized, designed nonlocal 
pseudopotentials are generated by the OPIUM package 
[33l l34| . The wave functions are expanded in a plane- 
wave basis with an energy cutoff of 680 eV. We initially 
consider the spin-orbit interaction based on a scalar- 
relativistic pseudopotential [Ejj and later, we will dis¬ 
cuss the effect of spin-orbit interactions, based on a fully- 
relativistic non-collinear scheme. 

The low-energy electronic structures of CU 3 NX are 
more or less similar for X = {Ni, Cu, Pd, An, Ag, Cd}, 
as reported in Ref. [21] ■ Here we present the results 
of Cu 3 NZn and Cu 3 NPd as representatives of transition 
metal-doped CU 3 N systems. Note that these are ex¬ 
treme cases where the transition metal atoms are maxi¬ 
mally doped [E]- In Cu 3 NZn the conduction and valence 
bands are mainly comprised of conduction A 2 U and va¬ 
lence Aig states near the Fermi energy. As shown in Fig. 
[^ these bands are inverted at the X points, forming two- 
dimensional (2D) Dirac points on the X-M and R-X 
lines (enclosed by red circles in the figure). These Dirac 
points signal the presence of a DLN enclosing X. Al¬ 
though there are more degenerate points near the Fermi 
level, and bands crossing the Fermi energy near the R 
point, we will simplify and here focus only on the bands 
near X. On the other hand, the conduction and valence 
bands of Cu 3 NPd are comprised of T 2 g and Ti„ states, 
which are inverted at the R point, forming the Dirac 
points on the R-X and M-R lines. These Dirac points 
are in fact parts of a DLN that encloses the R point, as 
shown below. 

The nodal lines of the conduction and valence bands 
in the 3D BZ are shown in Fig. [^ As mentioned above, 
DLNs appear near the X points in the Cu 3 NZn system. 
The cubic symmetry of the system dictates three DLNs 
encircling the three inequivalent X points A’' = nrja, 
where r = x,y, and z. Similarly, Cu 3 NPd also exhibits 
three DLNs due to the cubic symmetry, but since they 
appear enclosing the R point, they form in a gyroscope 
shape. In both systems, the DLNs are contained in three 
mirror-invariant planes at A®, X^ , and A^, due to the 
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FIG. 2. (color online) Electronic structures and Z 2 indices 
of (a) CuaNZn and (b) CusNPd. Bands are drawn along the 
high-symmetry lines of the BZ (inset). The Dirac points are 
indicated by red circles. Parity eigenvalues are illustrated at 
the eight parity-invariant points in the first octant of the BZ. 


corresponding mirror symmetries. We expect that break¬ 
ing the mirror symmetries should unlock the DLNs from 
the mirror planes, but that the DLNs will survive as they 
are protected by V and T. 

The appearance of DLNs agrees with the topologi¬ 
cal prediction of Z 2 invariants (vq] viV 2 Vs,)^ calculated 
from the parity analysis. In CuaNZn, parities at the 
eight time-reversal invariant momenta (L, ZX, 3M, R) 
give {vq\viV 2 V 3 ) = which dictates that there 

should be DLNs threading half the invariant plane at 
X'^ = -nr I a (r = x,y,z) an odd number of times. The 
three DLNs enclosing the X points fulfill this topologi¬ 
cal constraint (see the supplementary material for more 
details of this analysis). Similarly, in CusNPd we find 
that {vq; vih' 21 ' 3 ) = (1; 111), which is also in accordance 
with the formation of the three DLNs enclosing R. In 
this case, each invariant plane at X'^ is threaded three 
times by all three DLNs. 

A low-energy k ■ p Hamiltonian describing the conduc¬ 
tion A 2 U and valence Aig states, which form the DLNs 
in CuaNZn, captures the essential features of the DLNs. 
Near A’’’, symmetries dictate a two-band Hamiltonian 



(b) 



FIG. 3. (color online) Dirac line nodes in the Brillouin Zone 
(BZ). (a) and (b) GuaNZn, and (c) and (d) GuaNPd. The 
DLNs are illustrated by red curves in the 3D BZ [(a) and (c)] 
and on the 2D boundary plane of the BZ at fc = X^ [(b) and 

(d)]. 


A DLN forms at g,. = 0 and |q_Lp = Qq = —A/b±, when 
the bands are inverted (Ae < 0). The DFT results de¬ 
termine Ae ~ —0.4 eV. In CusNPd, unlike in CusNZn, 
there are conduction T 2 a and valence Tiu states, instead 
leading to a six-band Hamiltonian R. However, this can 
be decomposed into three copies of Hr with r = x,y, and 
2 ; and H — Hx®Hy®Hz, giving rise to three gyroscope¬ 
shaped DLNs. Therefore, the essential features of the 
DLNs should be the same between CusNZn and CuaNPd, 
aside from the former having a single DLN occurring in 
three inequivalent valleys of the BZ (X points) and the 
latter having three DLNs in a single valley (i? point). 

This model Hamiltonian also describes boundary 
modes. Consider a boundary perpendicular to f in which 
Ae varies between a negative (inverted) value and a large 
positive value. Fixing q± and considering the theory to 
linear order in qr —^ —idr, 

^z{q) = -wr^9r' + (Ae(r)-b6j_g^)r^-b(e-ba_Lgi)Ir. (5) 


Hr = {e+a±_\qi_\'^+arql)lT+vqrTy+{Ae+b±_\qi_\'^Hhrql)T'', 

(3) 

where q = k — X'", T represents the normal components 
to f, and the Pauli matrices describe the Aig 

and A 2 U states. The form of Hr is uniquely determined 
by inversion V = and time-reversal T = K (K being 
complex conjugation), together with the point group 
symmetries of X. It gives energy eigenvalues 

E±{q) = € + a±\q±\'^ + 

± {Ae + b±\q±\'^ + brq'^)'^ +v‘^q^. (4) 


For each k± this defines a Jackiw-Rebbi problem m- 
When Ae + b±k'j_ <0 there will be a boundary mode at 
the surface. In general the boundary band is not flat, but 
disperses for k± < fcp 

eo(fc_L) = e + < 0. (6) 

If a = 0 however, the surface band is flat. This reflects 
an additional chiral symmetry {H, t“} = 0 at this point. 
In this model, the value of a is related to the difference 
of the effective masses of the Aig and A 2 U bands. If the 
surface in the absence of inversion is electrically neutral. 
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FIG. 4. (color online) Two-dimensional surface electronic 
structure for CusNZn. (a) First octant of three-dimensional 
Brillouin zone (BZ) of CuaNZn projected onto the two- 
dimensional surface BZ of the (100) surface and (b) surface 
electronic band structure. The Dirac line nodes (DLNs) and 
the projected interior of DLNs are illustrated with the red 
and blue schemes, respectively. The slab bands are shown 
in black lines and surface states in the enclosed region are 
shown in blue lines. The shaded region represents bulk bands 
projected onto the BZ of the (100) surface along ■ 


then after inversion the surface will be neutral when the 
surface band is half-filled. This leads to a narrow surface 
band, where electron density = |Ae|/ 47 r 6 _L is con¬ 

trolled by the degree of band inversion. In the absence of 
screening from other bands, this surface band will tend 
to be pinned at the Fermi energy. 

To study the surface states in CuaNZn, we calculate 
the band structures of a slab geometry with 40 unit cells, 
exposing the (100) CU 2 N surfaces to vacuum. Our calcu¬ 
lation from first principles predicts that nearly-flat sur¬ 
face states emerge in the interiors of projected DLNs 
connecting the Dirac nodes, as shown in Fig. The 
slab band structure exhibits the weakly-dispersing sur¬ 
face states near F in the projected interior of the DLN. 
The topological surface states resulting from closed DLNs 
are half-filled and nearly flat, providing a unique venue 
for interesting strong-correlation and transport physics. 

The strong spin-orbit interaction can induce diverse 
topological phases in DLN semimetals, including topo¬ 
logical insulators, 3D Dirac semimetals 01351 , or even 
other DLN semimetals [13]. Analogously to graphene, 
spin-orbit interaction can gap out DLNs and drive the 
system to a topologically-insulating phase. The resultant 
topological insulator should have the same topological Z 2 
indexes as the DLN semimetal from which it originated. 
More interestingly, an additional crystalline symmetry 
may protect a part of the DLN in a symmetry-invariant 
region of the BZ, resulting in topological Dirac semimet¬ 
als or crystalline symmetry-protected DLNs with strong 
spin-orbit interactions. We have tested the effect of spin- 
orbit interaction in CusNPd using a fully-relativistic 
non-collinear scheme, and indeed found that C 4 sym¬ 
metry along the R-M line protects the Dirac point on 
the line, while the spin-orbit coupling otherwise opens a 
gap (with maximum size of ~ 62 meV on the R-X line), 
thus giving rise to a 3D Dirac semimetal phase in a strong 


spin-orbit interacting regime. Note that CusNPd is an 
extreme case where Pd is maximally doped, and thus 
the spin-orbit interactions due to Pd dd states are max¬ 
imized. The spin-orbit interaction can be controlled ei¬ 
ther by the Pd-doping concentration, or by doping other 
group-X transition-metal atoms, such as Ni, Pd, and Pt. 
We thus expect both the DLN semimetal and 3D Dirac 
semimetal phases should be accessible in the CU 3 N sys¬ 
tem. 

In summary, we have demonstrated that the combi¬ 
nation of inversion and time-reversal symmetries allows 
for the Z 2 classification of topological semimetals under 
vanishing spin-orbit interactions. The proposed topolog¬ 
ical semimetals are characterized by the presence of bulk 
DLNs and nearly-flat surface states, protected by inver¬ 
sion and time-reversal symmetries. Our first-principles 
calculations predict that the proposed topological phase 
can be observed in CU 3 N by doping with a class of non¬ 
magnetic transition metal atoms X, where X = {Ni, 
Cu, Pd, Ag, Cd}. The 2D surface states predicted for 
the DLN semimetal can hopefully be experimentally ob¬ 
served through, for example, ARPES in Cu 3 NXa;, using 
the doping concentration a; as a knob to control the sizes 
of the closed DLN and the enclosed surface band. Dop¬ 
ing with heavier atoms can also be used to potentially 
observe spin-orbit-induced topological phases. 

While this manuscript was in the final stages of prepa¬ 
ration we learned of recent work proposing DLN in Ca 3 P 2 
|39j . YK acknowledges support from NSF grant DMR- 
1120901. CLK acknowledges support from a Simons In¬ 
vestigator grant from the Simons Foundation. AMR ac¬ 
knowledge support from the DOE Office of Basic Energy 
Sciences, under grant number DE-EG02-07ER15920. 
Computational support is provided by the HPCMO of 
the U.S. DOD and the NERSC of the U.S. DOE. 


SUPPLEMENTARY MATERIAL FOR DIRAC 
LINE NODES IN INVERSION SYMMETRIC 
CRYSTALS 

Proof of Eqn. (1). 

Here we prove that the path integral of a Berry con¬ 
nection A(k) = -I X)„(Mn(k)|VkMn(k)) 




A-dk 


(SI) 


on a loop Cab = Cab — Cab that connects two parity- and 
time-reversal-invariant points F^ and F;, along the two 


time-reversal paths Cab and Cab (See Fig. SI I can be ob¬ 
tained by the parity eigenvalues ^„(Fa) = ±1 of the oc¬ 
cupied Bloch states at parity-invariant momenta Fq 


w(C'afe) = CaCfe; ■fa = ]^fn{rQ). 


(S2) 
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FIG. SI. Invariant loop Cab in the Brillouin zone (BZ). The invariant loop Cab = Cat — Cab connects two invariant points Fa 
and F;,, where Cab is a path from Fa to F;, along the time-reverse of Cab- 


It follows that 

uj{Cab) = 

= . (S3) 

Inversion symmetry (V) guarantees that |u„(—k)) = 
e*^"(*^)7^|M„(k)). It then follows that A(k) — A(—k) = 
VfcX]n/^n(k), so that 




(S4) 


Now consider P(k) = n„('“"(“k)|P|M„(k)) = 

/3n(k)^ Por a parity-invariant point F^ , P(Ti) = 
nn'fn(ri), where C«(ri) = (u„(rj)|P|'u„(ri)) is the par¬ 
ity eigenvalue. We thus prove Eqn. (S2). 


Z 2 topological invariants and Dirac line nodes of 
CugNZn 

The topological invariant = — 1 dictates that 

for any invariant surface Sabcd of the BZ, hosting four 
invariant momenta F^ (i = a,b,c,d), there will be an 
odd number of Dirac line nodes (DLNs) intersecting the 
half surface at k (and the other half at —k). Here we 
show that the DLNs that appear in CuaNZn satisfy this 
topological constraint. The cubic BZ of CusNZn has 
eight distinct invariant momenta (F, 3A, 3M, R), and the 
parity eigenvalues at the high-symmetry momenta are 
calculated as (^(F),^(A),^(M),^(i?)) = ( 1 ,— 1 , 1 , 1 ) as 
shown in Fig |S2| Similar to the Z 2 topological invariants 
of topological insulators m, 22 topological invariants 
(izq; ^ 11 ^ 2 ^ 3 ) can be defined in the DLN semimetals as 

(- 1 )"”= n (S5) 

rij — 0,1 

(S6) 

nj^i=0,l 


where Ci=(nin 2 « 3 ) parity eigenvalues at the eight in¬ 
variant momenta, Fj=(„j„ 2 n 3 ) = (nibi-|-n 2 b 2 -l-n 3 b 3 )/ 2 , 
with rij = 0 , 1 , and the primitive reciprocal lattice vectors 
bi. The Z 2 invariants in Cu 3 NZn are then obtained as 
(r'o; vii> 2 V 3 ) = ( 1 ; 111 ), which dictate that an odd number 
of DLNs will pierce half the nontrivial invariant surfaces, 


51467 , 52457 , and 53567 (See Fig. S2). 


First-principles calculations show that three distinct 
DLNs appear in Cu 3 NZn, contained in the boundary 
planes of the BZ as shown in Fig. |S3[ For convenience, 
we will refer to the DLN near the F^ as Li, where 
i = 1,2,3. The DLNs in CusNZn form in a manner 
that is consistent with the topological constraint, im¬ 
posed by {vq]ViV 2 V 3 ) = (1; 111). To show this, we first 
consider an L 3 that encloses F 3 in the BZ. As shown in 
Fig.[^ L 3 intersects three invariant surfaces, referred to 
as 53567 , 5 o 235, and 5oi36. The Z 2 topological invariants 
{vq',viV 2 V 3 ) = ( 1 ; 111 ) dictate that 53567 is nontrivial, 
and 5 o 235 and 5oi36 are trivial, so that the nontrivial 
53567 be pierced by an odd number of DLNs, while 
the trivial 5oi36 and 5o235 planes will be pierced by an 
even (including zero) number of DLNs. It is clear from 
Fig.j^that the trivial 5o235 and 5oi36 planes are pierced 
by two DLNs (L 3 and Li for 5oi36, and L 3 and L 2 for 
50235 ), and is thus consistent with ^a^b^c^d = 1, where 
(abed) = (0235), (0136). The nontrivial 53557 plane also 
satisfies the corresponding topological constraint. To 
show this, we construct a T-invariant loop C 35 that con¬ 
nects F 3 and F 5 avoiding the intersection with L 3 on 
the plane. The loop C 35 connecting F 3 and F 5 is bent 
down in the —k^ direction, while its time-reversed part¬ 
ner C 35 is bent up in the kz direction, as shown in Fig. 
S4 From the figure, it is clear that the invariant surface 
53567, containing is pierced once by L 3 in its half 
plane, thus satisfying the topological constraint imposed 
by ^ 3 ^ 5 ^ 6?7 = ~1- Therefore, the appearance of L 3 satis¬ 
fies the topological constraint. For the other inequivalent 
DLNs Li and L 2 , the cubic symmetry of Cu 3 NZn allows 
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FIG. S2. Six invariant surfaces Sabcd in CusNZn. (a) First octant of three-dimensional (3D) BZ. of CuaNZn are presented 
at parity-invariant momenta Fi, which determine the {vo', i'\U 2 Vz) = (1; 111) phase, (b) Shaded regions represent three trivial 
invariant planes (top panel) and three nontrivial invariant planes (bottom panel). 



FIG. S3. Dirac line nodes (DLNs) and invariant planes in the BZ of CuaNZn. Red circles represent the DLNs. The grey- 
shaded planes illustrate the invariant planes that the DLN at Fa point (Ls) intersects, one of which (S'ase?) is nontrivial and 
the other two (5o235 and Soiae) are trivial. The intersecting position is depicted by red scheme on the planes. 


us to apply the same argument for Li and L 2 to con¬ 
clude that all the DLNs in CusNZn obey the topologi¬ 
cal constraints characterized by Z 2 topological invariants 
(t'o; viV2Vz) = (1; 111). 


Band structures beyond CuaNZn and CuaNPd 

In this section we extend our discussion of the ma¬ 
terial realization of DLN semimetals. We demonstrate 
that DLNs can occur in CU 3 N by doping a more general 
class of non-magnetic Sd and Ad transition metals (TMs) 
atoms in the X, XJ, and XII groups of the periodic 
table, and that CuaNPda, can realize DLNs in a low- 
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FIG. S4. Invariant loop avoiding intersection with the DLN. The loop C35 connecting Fa and Fs is bent down in the —kz 
direction, while its time-reversed partner C35 is bent up in the fcz direction. The corresponding interior surface Sase? is threaded 
once by the DLN (illustrated by a red circle). 



FIG. S5. Band structures of GuaNX, with X={Vo, Ni, Zn}. (a) GuaN. (b) CuaNNi, and (c) GuaNGu. The Dirac nodes are 
indicated by red circles. 


doping concentration x (x < I). For this purpose, we 
first present the band structures of CU 3 NX, with X={Ni, 
Cu} in Fig. S5 These bands structures, including those 
of CusNZn and CuaNPd are more or less similar, when 
considering the position-shift of the Fermi energy due 
to the electron valence of the TM dopants. The Fermi 
energy positions one band higher in the group XII TM- 
doped case, comparing to that of the group X TM-doped 
cases, due to two more valence electrons in the group XII 
TMs than in the group X TMs. Comparing to the band 


structure of CU 3 N in Fig. S5 'a), it is clear that the group 
XII TMs provide two electrons per unit cell of CU 3 N, 
and thus move the preexisting DLNs near the Fermi en¬ 
ergy. This indicates that a high doping concentration x 
should be essential in the realization of the DLN near the 
Fermi energy in the group X//-doped systems. 

However, in the group X-doped systems, such as 


CusNPd, and Cu 3 NNi [see Fig. S5[b)], these high con¬ 
centrations {x ^ 1) are unnecessary. Under the doping of 
the group X TMs, the Fermi level remains in the same re¬ 
gion as where it was in CU 3 N, and instead doping drives 
the band-inversion. Therefore, DLNs start to appear 
near the Fermi energy in a lower concentration regime. 
To demonstrate this, we calculate the band structures of 
Cu 3 NPd„, with X = 0.0, 0.125, 0.25, 0.5, 0.75 and 1.0 by 
varying the number of Pd atoms in a 2x2x2 supercell 
of CU 3 N. The results, presented in Fig[^ show that the 


conduction and valence bands are already inverted even 
at a; = 0.125, and that the degree of the band inversion, 
as well as the size of the DLN, increase with increasing 
doping concentration x. 

Lastly, we consider the spin-orbit interaction in 
Cu 3 NPd. In Fig. |S7[ we present the band structures 
of Cu 3 NPd calculated with and without spin-orbit cou¬ 
pling (SOC). It is clear that the SOC induces a sizable 
gap in the DLNs up to ^ 62 meV along the R-X high- 
symmetry line in the BZ. However, the SOC cannot com¬ 
pletely gap out the entire DLN, as a single nodal point 
on the R-X line additionally protected by a rota¬ 
tional symmetry. The states forming the nodal point 
have distinct eigenvalues of the rotational symmetry op¬ 
eration, and thus retain the degeneracy even in the pres¬ 
ence of the SOC. This is in fact one of the mechanisms 
which stabilizes three-dimensional Dirac (3D) semimet¬ 
als O [38l So] • Therefore one might expect the 3D Dirac 
semimetal phase to be present in Cu 3 NPd. However, 
considering that the SOC mainly comes from the id or¬ 
bitals of Pd, and that the calculated SOC gap is esti¬ 
mated based on a somewhat extreme condition that Pd 
atoms are fully doped (one per unitcell of CU 3 N), it is 
more likely that the DLN semimetal should persist in a 
wide range of the Pd-doping concentration, and should 
be especially robust in the low-concentration regime. In 
order to strengthen (weaken) the SOC, one can substi- 





















































FIG. S6. Band structures of CusNPda:. (a) x = 0.0, (b) x = 0.125, (c) x = 0.25, (d) x = 0.5, (e) x = 0.75, (f) x = 1.0. The 
Dirac nodes are indicated by red circles. 



FIG. S7. Band structnres of CuaNPd with and without spin-orbit interaction. The bands in the rectangles in (a) are magnified 
in (b). 


tute the dopants from Ad TMs to 5d (3c?) TMs such as 
Pt (3d Ni). One can even dope with magnetic TMs to 
explore the effect of time-reversal-symmetry-breaking in 
the parent DLN semimetals. 
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